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■ Abstract 

We study the plasma-vacuum interface problem in relativistic magnetohydrodynamics for the 
case when the plasma density does not go to zero continuously, but jumps. In the vacuum region 
we consider the Maxwell equations for electric and magnetic fields. We show that a sufficiently 
large vacuum electric field can make the planar interface violently unstable. By using a suitable 
secondary symmetrization of the vacuum Maxwell equations, we find a sufficient condition that 
precludes violent instabilities. Under this condition we derive an energy a priori estimate in 
the anisotropic weighted Sobolev space for the variable coefficients linearized problem for 
nonplanar plasma-vacuum interfaces. 
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1 Introduction 

Consider the equations of relativistic magnetohydrodynamics (RMHD) for an ideal fluid (see, 

e.g., mm)- 

V a (pu a ) = 0, V a T afS = 0, V a (u a bP - uH a ) = 0, (1) 



where V a is the covariant derivative with respect to the Lorentzian metric g = diag (—1, 1, 1, 1) 
of the space-time with the components g a p Q p is the proper rest-mass density of the plasma, u a 
are components of the 4- velocity, 

T a ? = {ph + B 2 )u a uP + qg Ql3 - 6 a 6 /3 , 

h = 1 + e + (p/p) is the relativistic specific enthalpy, p is the pressure, e = e(p, S) is the specific 
internal energy, S is the specific entropy, B 2 = b a b a , b a are components of the magnetic field 
4- vector with respect to the plasma velocity, and q = p+ \B 2 is the total pressure. The 4- vectors 
satisfy the conditions u a u a = — 1 and u a b a = 0, and the speed of the light is equal to unity. 



We restrict ourselves to the case of special relativity, but in the future we hope to extend our results to general 
relativity (see discussion in Section [9|. 
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Let (x°,x) be inertial coordinates, with t = x° a time coordinate and x = (x , x 2 , x 3 ) space 
coordinates. Then 



u = —Uq = T, u 1 = Ui = Tvi, u := (u\, u 2 , u 3 ) = Tv, 

H; 



r 2 = l + |u| 2 , b° = -b = (u,H), b l = b i = ^ + (u,H)v l , 



b := (b x ,b 2 , b 3 ), B 2 = \b\ 2 - b 2 = 1^- + (v, H) 2 > 0, 

where T = (1 — |u| 2 ) -1 / 2 is the Lorentz factor, v = (1)1,1)2,1)3) is the plasma 3-velocity, H = 
(Hi, H2, H3) is the magnetic field 3- vector, the Latin indices run from 1 to 3, and by ( , ) we 
denote the scalar product. 

The RMHD equations ([!]) can be now rewritten in the form of the following system of 
conservation laws: 

d t ( P T) + div(pu) = 0, (2) 
d t (phFu + \H\ 2 v - (v,H)H) 

+ d\\ ((ph + B 2 )u®u-b®b) +Vq = 0, (3) 
d t (phF 2 + \H\ 2 -q) + div (phTu + \H\ 2 v - (v, H)H) = 0, (4) 
d t H- V x (vxH) = 0, (5) 

where dt = d/dt, V = (61,82,63), 8i = djdx 1 , etc. Moreover, the first equation in the Maxwell 
equations containing in JI]), 

divF = 0, (6) 

should be now considered as the divergence constraint on the initial data U(0,x) = Uq(x) for 
the unknown U = (p, u, H, S). 

Using flSJ) and the additional conservation law (entropy conservation) 

dt(pTS) + div (pSu) = (7) 

which holds on smooth solutions of system ([2"])- (p)|) , and following GODUNOV's symmetrization 
procedure we can symmetrize the conservation laws (|2"D-([5j) in terms of a vector of canonical 
variables Q = Q(U). This was done by Ruggeri k, Strumia [TU] and also by Anile & Pennisi 
[31 [2]. However, a concrete form of symmetric matrices was not found in |19t [3]. Moreover, if 
we deal with an initial-boundary value problem it is very inconvenient to work in terms of the 
vector Q. 

In principle, from Q we can return to the vector of primitive variables U keeping the symme- 
try property (see [S]). But, for RMHD finding a concrete form of symmetric matrices associated 
to the vector Q and returning from Q to U is connected with unimaginable (or even almost un- 
realizable in practice) calculations. To overcome this serious technical difficulty Freistuhler 
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& Trakhinin [8] symmetrized recently the RMHD system by using the Lorentz transform. 
Namely, taking into account ([6]) , for the fluid rest frame (v = 0) we can rewrite ([2]) , ([3|) , © , ([7|) 
in a nonconservative form which will be already a symmetric system. After that we should prop- 
erly apply the Lorentz transform to this system to get a corresponding symmetric system in the 
LAB-frame. Referring to [5] for details, here we just write down the resulting symmetrization 
of RMHD: 



A (U)d t U + J2M U ) d 3 U = °, 

3=1 



(8) 
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JT) — (v (8) ej + ej (8) u) 



r v J y r 

and I is the unit matrix. Clearly, A a are symmetric matrices. Note also that the last equation 
in © is just the nonconservative form dS/dt = of (J7J), with d/di = 9 t + (f, V). 

As is known [2], natural physical restrictions guaranteeing the hyperbolicity of the RMHD 
system do not depend on the magnetic field and coincide with corresponding ones in relativistic 
gas dynamics. In our case, by direct calculations one can show that the hyperbolicity condition 
Aq > holds provided that 

p > 0, p p > 0, < c 2 s < 1 (9) 

(of course, by default we also assume that \v\ < 1), where c s is the relativistic speed of sound, 
(? s = a? jh = Pp/h. The last inequality in ([9]) is the relativistic causality condition. 

Plasma-vacuum interface problems for the nonrelativistic MHD system usually appear in the 
mathematical modeling of plasma confinement by magnetic fields. Most of theoretical studies 
were devoted to finding stability criteria of equilibrium states, and the typical work in this 
direction is the classical paper [4J. According to our knowledge the first mathematical study 
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of the well-posedness of the full (non-stationary) plasma- vacuum model is the recent work |28j . 
This model consists in the MHD system in the plasma domain, the elliptic system V x % = 0, 
div% = (so-called pre-Maxwell dynamics) in the vacuum region, and corresponding interface 
boundary conditions [H[28], where by T~L we denote the vacuum magnetic field. 

The relativistic version of the classical plasma- vacuum interface problem [H [25J is of great 
interest in connection with various astrophysical applications such as, for example, neutron stars. 
Unlike the nonrelativistic case, we cannot now neglect the displacement current dtE and must 
consider the full Maxwell equations in vacuum: 

d t n + V x E = 0, d t E - V x U = 0, (10) 

where E is the electric field in the vacuum and we recall that the speed of the light is equal to 
unity. We do not include the equations 

div% = 0, div£ = (11) 

into the main system ()10p because they are just divergence constraints on the initial data 
V(0,x) = Vq{x) for the "vacuum" unknown V = (H,E). This is true for the Cauchy prob- 
lem, but below we will also prove this for our plasma-vacuum interface problem. 

Let us now discuss boundary conditions which we should pose on the interface (free bound- 
ary) between plasma and vacuum. Assume that this interface is a hypersurface 

= {x 1 = tp(t,x')}, with x' = (x 2 ,x 3 ). (12) 

Our assumption that interface f)12|) has the form of a graph is not a strong restriction for special 
relativity. Regarding the case of general relativity which we plan to consider in the future, 
the interface should be a compact codimension-1 surface. However, as for shock waves we can 
follow Majda's arguments |14] for extending the results to a compact plasma-vacuum interface 
(see also a related discussion in [27] for gas dynamics). On the other hand, even for general 
relativity, the above assumption is not still a strong restriction because, as in [7J, we can consider 
the interface not only locally in time (in the sense of a future local-in-time existence theorem), 
but also locally in space (for compactly supported initial data)@ 

Let ± (t) = jx 1 ^ ip(t, x')j are space-time domains occupied by the plasma and the vacuum 
respectively. That is, in the domain £l + (t) we consider system ©-([5]) (or (jSJ) governing the 
motion of an ideal relativistic plasma and in the domain 0~(t) we have the Maxwell equations 
(|10p in vacuum. The interface E(t) is a free boundary and moves with the velocity of plasma 
particles at the boundary: 

d t (f = v N on £(t), (13) 

where vn = {v,N) and iV = (1, — c^v?, —d^ip). 

Of course, locally we may suppose that the interface has the form of a graph. 
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The relativistic counterpart of the interface boundary condition [q] = g|s — ll^l^j = in 
[U [28] is the condition 

1-1/12 _|p|2 

- onS(t), (14) 



2 

Indeed, remembering that the plasma electric field E + = —v x H, the relativistic total pressure 
can be rewritten in the form 

'#| 2 {v,H) 2 \H\ 2 -\E+\ 2 



p + — ^r- H — = p + 



2r 2 2 ^ 2 

That is, the value (\H\ 2 — \E + \ 2 )/2 can be considered as the electromagnetic contribution (elec- 
tromagnetic pressure) to the relativistic total pressure in plasma, and in the right-hand side of 
equation (|14p we have the corresponding value in vacuum. Note also that the values \H \ 2 — \E + \ 2 
and \H\ 2 — \E\ 2 are Lorentz invariants, in particular, q = p + \\H' | 2 , where H' = H\ v= q is the 
magnetic field in the fluid rest frame (and E^_ = 0). 

As in the nonrelativistic case [U [28] , we have also the conditions 

H N = on S(t), (15) 

Hn = on E(t), (16) 

where T~L at = (^ , JV) and i^Ar = (H, N) . These conditions express the fact that the interface is a 
perfectly conducting free boundary. At the same time, as for current-vortex sheets [26J they are 
not real boundary conditions and should be considered as restrictions on the initial data (see 
Section [3] for the proof). 

At last, again as in the nonrelativistic case we have the boundary conditions 

Udt<p = N x E 

on S(t). However, the first of these boundary conditions is reduced to (|16|) . That is, we have 
finally the following two boundary conditions 

E 2 = n 3 d t ip - E^if, E 3 = -H 2 d t ip - E^ip on (17) 

for the vacuum electric and magnetic fields. It is worth noting that these boundary conditions 
are just jump conditions which can be deduced for the system b%H + V x = in Q^i), 
with E + = —v x H, E~ = E, H + = H, and H~ = T~L, i.e., for system ([5]) and the first system 
in (Tro|) , To get the final form (fT?]) we should take into account (fT3|) and (fT5|) , 

Thus, we obtain the free boundary vale problem for system (jSJ) in Q + (t) and system (|1U|) in 
f2~(i) with the boundary conditions (fT3|). ([14^) . (fT7|) on £(t) and the initial data 

U(0,x) = U (x), x£n + (0), V(0,x) = V Q (x), iefi-(0), 

(18) 

tp(0,x') = tp (x'), xeE 2 . 
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Moreover, we will prove that ©, (fTT|h (fT5|) and (fT6j) are restrictions on the initial data (fTBj) . i.e., 
if they are satisfied at the first moment t = 0, then they hold for all t > 0. System (|10p is always 
hyperbolic and we assume that the hyperbolicity condition ([9]) is satisfied up to the boundary 
of the domain £l + (t). That is, as in the fluid- vacuum problem in |12[ 127] . we consider the case 
when the density does not go to zero continuously, but jumps. However, unlike the problem in 
[12\ 127] . we do not have the boundary condition p\-£ = and, therefore, the condition p|s > 
does not contradict, for example, the equation of state of a polytropic gas. We may suppose 
that the jump [p] = p\s is very small that corresponds to a small but non-zero pressure on the 
interface. 

Our final goal is to find (at least sufficient) conditions on the initial data (|18p providing the 
local- in-time existence and uniqueness of a solution (U, V, ip) of problem (JSj) , (fTUj) . (fT3j) . JHI), (fT7|) . 
(|18p in Sobolev spaces. The first step towards the proof of such an existence theorem is a careful 
study of a corresponding linearized problem. The main result of the present paper is finding a 
sufficient stability condition which gives a (basic) energy a priori estimate. Roughly speaking, 
we show that if the unperturbed vacuum electric field is small enough and the unperturbed 
plasma and vacuum magnetic fields are nowhere parallel to each other on the interface, then 
the interface is stable. It is noteworthy that this condition could be essentially specified by a 
numerical analysis of the positive definiteness of some matrix of order 42 (see (|126p in Section 
[6]) . The crucial role in deriving a basic a priori estimate for the linearized problem is played by a 
secondary symmetrization of the Maxwell equations (|10p (see Section [2]) . Using the terminology 
of [25], we construct a dissipative 1-symmetrizer for the linearized problem. This allows us to 
find the mentioned stability condition. Moreover, we show that the corresponding constant 
coefficients linearized problem for a planar interface can be violently unstable for a sufficiently 
large vacuum electric field. 

The main difficulty in the study of the relativistic plasma- vacuum problem is the same as for 
compressible current- vortex sheets \24\ I26j for the MHD system, i.e., we cannot test the Kreiss- 
Lopatinski condition |10|, 1 14j analytically. On the other hand, since, as for current- vortex sheets, 
the number of dimensionless parameters for the constant coefficients linearized problem is big, 
a complete numerical test of the Kreiss-Lopatinski condition seems unrealizable in practice. We 
overcome this principal difficulty by using the energy method and constructing the mentioned 
dissipative 1-symmetrizer. 

Another principal difficulties are again the same as for current-vortex sheets. That is, the 
Kreiss-Lopatinski condition can be satisfied only in a weak sense and the plasma-vacuum inter- 
face is a characteristic free boundary. Therefore, we have a loss of derivatives in our a priori 
estimates and a natural loss of control on derivatives in the normal direction that cannot be 
compensated in RMHD (unlike the situation in gas dynamics [3(27]). Therefore, in our nonlin- 
ear analysis (this work is now in progress) we have to use a suitable Nash-Moser-type iteration 
scheme (as, for example, in [61 [261 12Z]) and the natural functional setting is provided by the 
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anisotropic weighted Sobolev spaces H™ (see [29, 22J and Section 2] for their definition). In this 
paper, we prove our basic a priori estimate in the space H\ for the "plasma" unknown U, and 
for the "vacuum" unknown V we have a full regularity in the normal direction and can work in 
the usual Sobolev space H . 

The rest of the paper is organized as follows. In Section [2] we write down the mentioned 
secondary symmetrization of the vacuum Maxwell equations. In Section [3] we reduce our free 
boundary problem to an initial-boundary value problem in a fixed domain and discuss properties 
of the reduced problem. In Section SI we obtain the linearized problem and formulate the main 
result for it (Theorem 14. ip . In Section [5] we properly reduce the linearized problem to that with 
homogeneous boundary conditions and homogeneous Maxwell equations. Section \6\ is devoted 
to the derivation of the basic a priori estimate for the reduced linearized problem in the case of 
constant coefficients. In Section [TJ by considering particular cases of the unperturbed constant 
solution, we prove that the constant coefficients linearized problem can be ill-posed. In Section 
[8j we extend the result of Section[6]to the case of variable coefficients. At last, Section [9] contains 
concluding remarks. 



2 Secondary symmetrization of the vacuum Maxwell equations 



The Maxwell equations (jlOp form the symmetric hyperbolic system 



8 t V + Y, B AV = 0, 
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At the same time, they imply the wave equation for % (and E) which can be, in turn, sym- 
metrized, i.e., it can be rewritten as a symmetric hyperbolic system for derivatives of T~L (and E). 
There are different ways to symmetrize the wave equation, and symmetrizations may depend on 
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arbitrary constants or functions (see, e.g., [3 [25]). However, for our present goals it is better 
to use a counterpart of one of these symmetrizations for the original system (|19|) . Since the 
original system was already symmetric, it is natural to call its new symmetric form secondary 
symmetrization. 

Consider (|19p in the whole space M 3 . Clearly, 



d 
dt 



\V\ 2 dx = 0, 



but, using constraints (llip . by simple manipulations from (llOj) we can deduce the following 
additional conserved integrals: 

d 



d 
dt 



Hence, we have 



f (H 2 E 3 -n 3 E 2 )dx = 0, 
dt 7m3 

I (U 3 E 1 -H 1 E 3 )dx = 0, -[ (H 1 E 2 -U 2 E 1 )dx = 0. 
Jm.3 dt J R 3 

4 / {\V\ 2 + vi(H 2 E 3 - U 3 E 2 ) + v 2 {H 3 E 1 - H X E 3 ) + u 3 (HiE 2 - n 2 E±)} dx = 0, (20) 
dt J R 3 

where V{ are yet arbitrary constants, but below they will be arbitrary functions satisfying a 
certain hyperbolicity condition. 

The energy identity ([20|) also reads 



dt J R 3 



dx = 0, 



where 
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We now assume that vi in Bq are arbitrary functions Vi(t,x). Then, taking into account the 

divergence constraints (fTTj) . it follows from (fT§|) that 

3 3 
B d t V + BoBjdjV + i?idiv U + i? 2 div E = B d t V + ^ BjdjV = 0, 



(21) 
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The symmetric system (|2~Tj) is hyperbolic if So > 0, i.e., 

\v\ < 1, (22) 

with the vector- function v = (ui, ^3). Following [26], we can call the set S = (So, R±, R2) the 
secondary generalized Friedrichs symmetrizer. 

We will use (|2ip for obtaining a dissipative energy integral for the linearized plasma- vacuum 
problem, and for further arguments we will need the representations 

Bj = BqBj + fCj (23) 

following from the fact that 

3 

Rid\vU + R 2 divE = ^K.jdjV, JCj =I 2 ®(v<g) ej ), 
where I 2 is the unit matrix of order 2. 

3 Reduced problem in a fixed domain and its properties 

We now return to our free boundary value problem (jHJ), (fl3"]l . (JHJ), dTTJ) (JT9J) and reduce it to 
an initial-boundary value problem in the fixed space domain R% = {x 1 > 0, x' £ M 2 }. In 
[26\ [271 128j . to avoid assumptions about compact support of the initial data in the existence 
theorem and work globally in M$_ we used a cut-off function xi xl ) £ Co°(^) m ^ ne cnan g e of 
independent variables. Since in the future we plan to extend our results to the general relativistic 
version of the plasma- vacuum interface problem, here we prefer to apply a simplest straightening 
of the interface S (without using the cut-off x)- That is, in our future local-in-time existence 
theorem we will consider compactly supported initial data as, for example, in [3(7]. 

We make the simple change of variables x 1 = x 1 — (p(t, x') for the RMHD system (JSj) and 
x 1 = —x 1 + <f(t,x') for the Maxwell equations (|1U|) . In other words, the unknowns U and V 
being smooth in Q^ft) are replaced by the vector- functions 

U(t, x) := U(t, x 1 + <p(t, x'), x'), V(t, x) := V(t, -x 1 + tp(t, x'),x'), 

which are smooth in the half-space Dropping for convenience tildes, we reduce ([8]), f)13|) . 
(HH), (fT7|) ~ ^9|) to the initial boundary value problem 

F(U,<p) = in [0,T] x R%, (24) 
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Y(V,(p) = in[0,r]xR^_, 
T, V,(p) = on [0, T] x {x 1 = 0} x R 2 , 
(U,V)\ t=0 = (U ,Vo) mRl, <p\t=a = <Po m M 2 , 
where F(U,tp) = L(U,<p)U, V(V,ip) = M(cp)V, 

L(U, if) = A (U)d t + M{U, <p)d! + A 2 {U)8 2 + A 3 (U)d 3 , 

MU, y>) = A 1 (U) - A (U)d t <p - A 2 (U)d 2 ip - A 3 (U)d 3 ip, 
M{y) =Id t - B x {(p)dx + B 2 d 2 + B 3 8 3 , B x {<p) = B 1 -Id t <p- B 2 d 2 ^ - B 3 d 3 <p, 
( v N - d t ip \ 



(25) 
(26) 
(27) 



l(U,V,<p) 



Regarding condition (|15 



)(\n\ 2 - \e\ 2 ) 
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\ E T3 + n 2 d t if j 



, E 7i =E 1 8 i <p + E i , i = 2,3. 



HnL^o = 0) ( 28 ) 
and the divergence constraint Q, which in the straightened variables takes the form 



div h = 0, 

with h = (J3jv, H 2 , H 3 ), referring to |26| for the proof, we have the following proposition. 



(29) 



Proposition 3.1 Let the initial data ((27]) satisfy fl28|) and P?]) . If problem (f2^|) - (|27|) has a 

solution (U,V,<p), then this solution satisfies (|'28|) and (|29p for all t £ [0,T]. 

In the proof of an analogous fact (see Proposition (|3.2p below) for the condition 

H N \ x i =0 = (30) 



(31) 



and equations (jlip written in the straightened variables, 



div \] = 0, div e = 0, 



where f) = (Hn, W-2, W3), e = (En,E 2 ,E 3 ), E n = (E,N), and div a = -dia\ + d 2 a 2 + d 3 a 3 for 
any vector a = (01,02,03), we should distinguish between two cases. In the first case 

dtip < (32) 

and the strict inequality in (|32l) corresponds to the plasma expansion. In the opposite case 

d t ip > (33) 
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we have shrinkage of the vacuum region. For the most general mixed case when for some parts of 
the interface the plasma expands into the vacuum and for other parts the vacuum expands into 
the plasma we below show that the boundary x 1 = for the Maxwell system (|25p is of variable 
multiplicity. We do not consider this very difficult case. At the same time, since, as already 
noted in Sect. 1, in our future local- in-time existence theorem we will consider the interface 
locally in space (for compactly supported initial data), it is reasonable to assume that one of 
the conditions in (|32|) and f|33|) is satisfied for the whole interface. 

We now show that the boundary matrix —B\{ip) for system (|25p has different signatures for 
cases ([32]) and ([33]) . We can directly calculate its eigenvalues: 

Ai,2 = d t ip + y/l + (d 2 tp) 2 + (d 3 ip) 2 , A 3>4 = d t ip - ^1 + {d 2 p) 2 + {d^) 2 , A 5i6 = d t cp. 

Since the interface speed should be less than the speed of the light, i.e., 

_ M _ < 1 

the eigenvalues Ai 5 2 > and < for case ([32]) and Ai^ < and A3_4 > for case ([33]) . Note 
that we even have the more restrictive inequality 

\dt<p\ < 1 (34) 

because the first boundary condition in ([26]) for the particular case V2\ x i=q = V3\ x i =0 = reads 
dtf = vi\ x i =0 . 

Let us count the number of incoming characteristics for systems (|24|) and (|25p . The simplest 
way to do this for the RMHD system (|24|) is to write down the quadratic form with the boundary 
matrix Ai(U, ip) on the boundary x 1 = computed on a given basic state (U, (p) ("unperturbed 
flow") about which we will linearize system (]24p . Moreover, later on we will also need this 
quadratic form for deriving a priori estimates for the linearized problem by the energy method. 
We assume that the basic state satisfies (]28p and the first boundary condition in (|26p . In fact, 
because of this assumption the boundary x 1 = is characteristic for system (j24|) and we now 
prove this. 

It is convenient to use the following representations for the symmetric matrices Aj in ([8]): 

Aj = VjA + Gj, 

( ej-v jV T ^ 

ej - v jV Qj Mj T 
V J\fj J 



with 

Gj = 
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(v H) B 2 
+ ' (.ff <8> ej + ej <g> H) - — (v <g> ej + ej <g> v) 



Then, taking into account the first boundary condition in (|26p . we have 

.4 := MU, f>)\ x i =0 = Gt(U)\ x i= - G 2 (U)\ x i =0 d 2 v - G 3 (U)\ x i =o d 3 0. 
Using ([28]) , after some algebra we obtain 

(AU,U) = 2Tqv N \ x i =0 , (35) 

where 

v N = {v,N), N=(l,-d 2 <p,-d 3 £), q=p+^SB 2 , 
v is the perturbation of u/F, and 5B 2 is the perturbation of the electromagnetic pressure B 2 : 
v = i(u- (v,u)v), 5B 2 = ^((b,H) + (v,H)(H,u)-B 2 (v,u)y (36) 

Here all the "hat" values are calculated on the basic state, e.g., v = u/F, b = H /F + (u,H)v, 
etc. 

Let W = (q,Fv N ,u 2 ,u 3 ,H,S) and U = JW . Clearly, det J / 0. It follows from ([35]) that 

/ 1 ••• \ 



(AU, U) = {J T AJ W, W) = (S 12 W, W), S 
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1 




(37) 



y ••• j 

Hence, the matrix A has one positive ("incoming") and one negative eigenvalue, and other 
eigenvalues are zeros. 

In case ([32]) system ([25]) has two positive eigenvalues, and the whole system (|24p . ([25]) for 
U and V has three incoming characteristics. In this case the number of boundary conditions in 
(|26|) is correct because one of them is needed for determining the function (p(t,x'). If, however, 
the vacuum region expands, i.e., condition (|33[) holds, than the correct number of boundary 
conditions is six, but we have only four in (|26p . It means that problem (|24p -(|27 p is formally un- 
derdetermined and we miss two boundary conditions. However, if as these two missing boundary 
conditions we take 

div~ fj |j.i = o = and div~ e| x i =0 = 0, (38) 

then, as for case (j32p . we can show that equations ([31 p are just divergence constraints on the 
initial data (|27p . Namely, the following proposition takes place. 
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Proposition 3.2 Let the initial data ([27} satisfy ([H and (|3Tj> . If problem (lMj) -(l27 |) /ias a 
solution (U,V,(p) with property ()32j) . i/ien i/iis solution satisfies (|30p and f)31[) /or a// i £ [0,T]. 
If problem (|24 p - (|27p to£/i £/ie too additional boundary conditions (|38p /ias a solution (U,V,ip) 
with property (|33p . i/ien i/iis solution again satisfies (|30p and (|3ip /or a// £ E [0,T]. 

Proof. System (|25p can be rewritten as follows 

$t) + V-x e + %>&f) + 5(^ = 0, (39) 
3te- V" xi5 + %>0j.e + .B(<p)e = O, (40) 

where 

15 = (n!,n T2 ,n^), <e = (e 1 ,e T2 ,e T3 ), 

( dtc^y? <9f<9 3 ^ ^ 

n n =n 1 d l ip + n l , 1 = 2,3, b = b{v)= 

\ / 

and the operator V~x differs from the usual curl operator only by the fact that V - = {—d\, d 2 , 
83). Applying the operator div~ to ([39]) and (|4"0"|) . we get the equation 

<9 t F + dtipdxF = in [0, T] x 

for i 7 = div - f) and .F = div - e. In case (|32p this equation has no incoming characteristics 
and, therefore, it does not need boundary conditions. Hence, if F\t=o = 0, then F = for all 
t G [0,T]. This yields ([3T]) . In case ([33]) . in view of the boundary conditions ([38]) . we have 



F\ 



x l =0 



0. Then, we again deduce ([3"T]) if F| t= o = 0. 



The first equation in (|39p reads 

5 t ^jv + 9i^tv + d 2 E T3 - d 3 E T2 + H 2 d t d 2 ^ + H 3 d t d 3 <p = 0. 

Considering it on the boundary x 1 = and using the last two boundary conditions in (|26|) . we 
obtain 

dtH-N — <9tv?div~ f) = at x 1 = 0. 

In view of ([3T]) . this implies e?t%jv|a;i=o = 0- That is, we get ([30]) for all t G [0, T] if "HjvUi=o = 
for t = 0. □ 

Below without loss of generality we consider the case (|32p . The proof of the basic a priori 
estimate for the linearized problem in Sect. [6] and [8] for the opposite case (|33p stays the same 
as for ([32]) . However, since the additional boundary conditions ([38]) are not quite standard, the 
preparatory arguments in Sect. should be modified and the future proof of the existence of 
solutions will be diferent from that for case (|32|) . Anyway, in this paper we restrict ourselves to 
case (|32p . Moreover, it is natural to assume that the "stable" counterpart 



d t (p < 



(41) 
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of condition f)32[) holds. In particular, in our future local-in-time existence theorem we will 
consider initial data satisfying condition (|4ip . We note also that under assumption (j41[) the 
boundary x 1 = is noncharacteristic for system ()25p Pl 

4 Linearized problem and main result 

Let 

(U(t,x),V(t,x),tp(t,x')) (42) 
be a given sufficiently smooth vector-function with U = (p, u,H,S), V = (H, E), and 

\\(U,V)\\wl(n T ) + \\di(U,V)\\wUn T ) + WWw&VOr) < K, (43) 
where K > is a constant, 

Q T ■- (_oo, T] x R 3 + , dtt T ■= (-oo, T] x {x 1 = 0} x M 2 , 

and below all the "hat" values will be related to the basic state (|34p . e.g., v = u/T, T 2 = 1 + |n| 2 , 
a 2 = l/pp(p,S), etc. 

We assume that the basic state (|4"2"j) satisfies the hyperbolicity condition Q in fL?, 

p(p,S)>0, Pp (p,S)>0, 0<c 2 <l, (44) 

together with such natural assumptions as \v\ < 1, inequality (jM)) for <p, etc. We also assume 
that P2|) satisfies the boundary conditions in ([26]) except the second one on 

vn\x x =o = x : E T2 \ x i =0 = xii 3 \ x i =0 , E T3 \ x i =0 = -X^ 2 |x 1 =0> ( 45 ) 

and the interior equations for H and % in containing in (|24p and (|25p: 

a t # + (u;,V)#-(h,V)i) + -H'divu; = 0, (46) 

Sfcjj + xdijj + V-xe + BC^^O, (47) 

where x = dt<p, w = (vn — ^2, ^3)- At last, we assume that the basic state satisfies condition 
dSJ). More precisely, let 

x < -5 < (48) 

on dQ-T with a fixed constant 5. Recall also that |x| < 1, cf. ()34p . 

Later on, for the linearized problem we will need equations associated to the nonlinear 
constraints (|28j) . (i29j) . (|30|) . and (|3T]l . However, to deduce them it is not enough that these 
constraints are satisfied by the basic state (|42p . For the plasma magnetic field, as in [26], we 

3 Even if dt<p = 0, the boundary is characteristic for (|25[) only formally because, in view of 1)31 p. we have the 
full control on derivatives of V in the normal direction. 
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need (JM|), and for the same reason for the vacuum magnetic field we have to use ([37]) . If we 
assume that the basic state satisfies 

^Jv|a;i=0 = 0' ^ArUi=0 = °> (49) 

div h = 0, div~ jj = (50) 

for t = 0, then it follows from the proofs of Propositions 13. ll (see |26| ) and 13.21 that conditions 
(|49p are true on 8VLt and equations (|50p hold in fi^. Thus, without loss of generality for the 
basic state we also require the fulfillment of (|49p on 8Qt and (|50|) in f^. 

Conditions ([49 P and the first condition in ([45 P endow the boundary matrix 

diag (2i(f7, 0), Ui=o 

for the linearized problem with the same structure as for the boundary matrix at x 1 = for the 
nonlinear problem ([24J) - ([27p . Concerning the second and third conditions in (|45|) . we will need 
them for deducing a linearized version of (|30p for solutions of the linearized problem. 

Remark 4.1 Assumptions (|44[) — (|48j) are nonlinear constraints on the basic state. In the forth- 
coming nonlinear analysis we plan to use the Nash-Moser method. As in [5^ 126} 127], the Nash- 
Moser procedure will be not completely standard. Namely, at each nth Nash-Moser iteration step 
we will have to construct an intermediate state (£^1+1/21^+1/2)^+1/2) satisfying constraints 
(|44|) -([48 p (and ([49 1) . (|5U|) for t = 0). For the "plasma" part of the problem such an intermediate 
state can be constructed in exactly the same manner as in [26J and for the "vacuum" component 
V n+ i/ 2 the construction is also not difficult. We omit corresponding arguments and postpone 
them to the nonlinear analysis. 

The linearized equations for ()24|) — (|26 [) read: 

P\U, 0)(SU, Sip) := ^-F(U e , <p £ )\ £=0 = fi in U T , 
de 

Y(V, $)(6V, Sip) := ^V(K, c^)| £=0 = fu in n T , 

M'(U,V,<p){5U,SV,5<p) :=-^B(U £ ,V £ ,ip £ )\ £=0 = g on 8Q T , 

de 

where U £ = U + eSU, V £ = V + eSV , and ip £ = (p + sSip. Here we introduce the source terms 
/ = (/i> /11) = (/i) • • • j fu) an d 9 = (#1 1 • • • j 9i) t° make the interior equations and the boundary 
conditions inhomogeneous. We easily compute the exact form of the linearized equations (below 
we drop 5): 

¥'(U, <p)(U, <p) = L(U, p)U + C(U, p)U - {L(U, 0)ip}d x U = f h 
V(V, <p)(V, p) = M(p)V + {M{$)<p}dxV = fu, 
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V(U, V,<p)(U,V,<p) 



( vn - d t (f - v 2 d 2 p - v 3 d 3 p \ 
q-(H,H) + (E,E) 
E T2 - H 3 d t <p - xU 3 + E x d 2 p 
\ E T3 + U 2 d t p + xU 2 + Eid 3 p ) 



9, 



where 



q = p + 5B 2 , v N = vi - v 2 d 2 p - v 3 d 3 p, E T . = E\dj(p + Ej, 



v and 5B 2 are written down in (J36|) . and the matrix C(U, tp) is determined as follows: 

C(U, tp)Y = (Y, X7 y A (U))d t U + (Y, V y A x {U, 0))d 1 U 

+ (Y, V y A 2 {U))d 2 U + (Y, V y A 3 (U))d 3 U, 



(Y,V y A(U)):=± yi 



Y=U 



, Y = (2/1, • • • ,Vs)- 



The differential operators W(U, (p) and V'(V, <p) are first-order operators in (p that gives some 
trouble in obtaining a priori estimates. As in [6j 124 1 f2G | I27j. following Alinhac pQ, we introduce 
the "good unknowns" 

U :=U -<pdiU, V :=V + pd{V. (51) 

Omitting detailed calculations, we rewrite our linearized interior equations in terms of the new 
unknowns (|5ip : 

L(U,p)U + C(U,p)U + pd 1 {F(U,p)} = f h (52) 

M(0)V-cpd 1 {W(V,0)} = f lI , (53) 

As in [TJ |6l [26j 127] . we now drop the zero-order terms in tp in (|52p and (|53|l R and consider 
so-called effective linear operators 

¥' e (U,p)U := L{U,p)U + C(U,p)U and Y' e {p)V := M (p)V . 

We write down the final form of our linearized problem for (U, V, p): 

L(U,p)U + C(U,p)U = fi mQ T , 
M{p)V = f u in Q T , 

( v N - d t p - v 2 d 2 p> - v 3 d 3 p + ip dii N \ 
q-(H,H) + (E,E) + [d iq }p 
E T2 - d t (n 3 p) - xii 3 + d^p) 

\ e T3 + d t (n 2 p) + ku 2 + ds(Ex<p) J 



g on 8Ut, 



(54) 
(55) 



(56) 



4 In the future nonlinear analysis the dropped terms in (|52p and (| 53p should be considered as error terms at 
each Nash-Moser iteration step. 



16 



(U,V,cp) = Q fort<0, (57) 

where 

vn = vi- v 2 d 2 (p - vzdzy, E T , = Etdjip + Ej, 

[diq] = (3ig)U =0 - (£,di£)U =0 + (E,diE)\ xl=0 , 

and v and q are determined similarly to v and q. We used the second and third equations in 
(|47p taken at x 1 = while writing down the last two boundary conditions in (|56p . We assume 
that / and g vanish in the past and consider the case of zero initial data, which is the usual 
assumption. We postpone the case of nonzero initial data to the nonlinear analysis (construction 
of a so-called approximate solution; see, e.g., |26j). 

Before formulating the main result for problem (|54 [ ) — f)5T[) we give the definition of the 
anisotropic weighted Sobolev spaces H™. The functional space H™ is defined as follows (see 
[291 (13 122] and references therein): 

H™(R S + ) := G L 2 {Ml ) \ d%d\u G L 2 (R 3 + ) if \a\ + 2k < m } , 

where m G N, d° = (adi^d^d^ 3 , and a(x 1 ) G C°°(IR+) is a monotone increasing function 
such that ^{x 1 ) = x 1 in a neighborhood of the origin and c(x x ) = 1 for x 1 large enough. The 
space H™(M.+) is normed by 



mi:,,.,- Y, W d * d l u \\l, 



\a\+2k<m 



We also define the space 



fc=o 



HT^t) = H H k ((-oo,T],H™- k (M.l )) 

equipped with the norm 

/i "- 
\\\u(t)\\\ 2 m ^dt, where = 
-OO „■_ n 



T 

2 

m-j',*- 

i=o 



Within this paper we use the space H™(Qt) mainly for m = 1. Clearly, the norm for HI(Qt) 
reads 



/ (u 2 + (^u) 2 + (d9iu) 2 + {d 2 uf + (d 3 u) 2 ) dtdx. 
We are now in a position to state the main result of this paper. 

Theorem 4.1 Let the basic state ([4*2]) satisfies assumptions (pl3|) - ([50|) . Zei a/so 

\H 2 ii 3 - H 3 n 2 \ lxl=Q >e>0, (58) 
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where e is a fixed constant. Then there exists a positive constant \i* such that for all (/, g) £ 
H*(£It) x -ff 3 (<9$7-r) which vanish in the past and for all fi < fj,*, where 

fi = \Ei+v 2 H 3 -v 3 H 2 \ ]xl=Q , (59) 

the a priori estimate 

[U]i,*,t + ||^||ifi(n T ) + IMItf^nr) ^ c {[fih,*,T + ll/ii||^3(n T ) + \\g\\ H 3(dn T )} ( 60 ) 

holds for problem (j54j) (|57p . where C = C(K,T) > is a constant independent of the data 
(/,<?)• 

5 Reduction to homogeneous boundary conditions and homo- 
geneous Maxwell equations 

Before deriving the basic a priori estimate (|60p we discuss some useful properties of the linearized 
problem (|54 p -(|57 p and reduce it to a problem with fu = and g = 0. First of all, as for current- 
vortex sheets [26J, we can prove the following proposition (see Appendix A in [26] for its proof 
for a technically more difficult case when the cut-off function x( xl ) is used in the straightening 
of the front £). 

Proposition 5.1 Let the basic state (|42p satisfies assumptions ()43|) — f|50 p . Then solutions of 
problem ([54"]) - (j5T]) satisfy 

div h = /15 in $7^, (61) 

H N - H 2 d 2 ip - H 3 d 3 (p + <p d\H N = g 5 on dVtr- (62) 

Here h = (Hn,H 2 ,H 3 ), Hn = H\ — H 2 d 2 (p — H 3 d 3 0, the new source terms /15 = fi^{t,x) and 
95 = 95(ti x ')> which vanish in the past, are determined by the source terms fe, fy, ft, g\ and 
the basic state as solutions to the linear inhomogeneous equations 

d t fi5 + (w, V/15) + /i 5 div w = T in fi T , (63) 

dt95 + v 2 d 2 g 5 + v 3 d395 + (d 2 v 2 + d 3 v 3 ) g 5 = Q on dtl T , (64) 

where 

T = div f h , f h = (f N , f 6 , f 7 ), f N = f 5 - f 6 d 2 ip - fid 3 <p, 
Q = {f N + d 2 (H 2 gi) + d 3 (H m ) } \ xl=Q . 

For deriving estimate (j60|) we will also need a "linear" version of Proposition 13.21 for the 
"vacuum" part of problem ([54"]) -(|57 p . 
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Proposition 5.2 Let the basic state (|4"2"]1 satisfies assumptions ([4"3]) - ([50]) . Then solutions of 
problem ([54"]) - ([57]) satisfy 

div~ f) = /i6, div~ e = /17 in T , (65) 

— H 2 d 2 <p — H 3 d 3 ip + (pd\7iN = 96 on dflx- (66) 

#ere f) = (H N ,H 2 ,H 3 ), H N = Hi - H 2 d 2 - H 3 d 3 (p, e = (E N ,E 2 ,E 3 ), E N = E 1 - E 2 8 2 - 
E 3 d 3 (p, the new source terms fi% = fi§{t,x), fn = fn(t,x) and g§ = g%(t,x'), which vanish in 
the past, are determined by the source terms /g,... , /14, 53, #4 and 6asic state as solutions 
to the linear inhomogeneous equations 

9tfm+k + xdihe+k = Fi+k (A; = 0,1) in tt T (67) 
dtge = Ql on dtt T , (68) 

where 

Fi = div~ f H , f n = (f}f, /10, /n), fh = h~ hod 2 - fnd 3 <p, 
T 2 = div~ /b, /b = (/l^ /13, /i 4 ), = / 12 - /i 3 <9 2 ^ - /u^, 

£l = (/A - #294 + $333 + /l6)U=0- 



Proof. Introducing 

h = {Ui,U T2 iU T3 ) and £ = (E 1 ,E T2 ,E T3 ) 
with H T . = %idi(p + Hi, we rewrite system ([56]) as (cf. ([39]) . PU]) ) 

fltfj + V- x 6 + x0i& + fl(0)& = / w , (69) 
d t e- V" xSS + xdii + B(0)t = f E . (70) 

Applying the operator div~ to ([69]) and ([TO]) , we obtain equations ([67]) for /ig = div~ rj and 
fn = div~ e. Considering the first equation in (|69p at x 1 = and using the last two boundary 
conditions in ([ST]) and the first equation in ([6"5"]h we get ([6"8"]) for 

g 6 =H N -d 2 (H 2 <p)-d 3 (H 3 <p). (71) 

Taking into account the second condition in (]50p on 30^, we obtain that the right-hand side of 
(|71|) coincides with the left-hand side of equation ([BT?]) . □ 

It follows from the first condition in (]45p and assumption (]48|) that the interior equations 
(|64p and ([67|) do not need boundary conditions. Therefore, using standard arguments of the 
energy method, from (|64p and (|67p we deduce the estimates 

II/i5||l 2 (q t ) < CII/iIIhH^t) - C[/i]v,t, (72) 
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IK/l6,/l7)||H 2 (f7 T ) + \\(fw, fn)\x 1 =o\\H 2 (dQ T ) - C \\M\h 3 {Q t )i (73) 

Here an later on C is a constant that can change from line to line, and sometimes we show the 
dependence of C from another constants. In particular, in (|72p the constant C depends on K 
and T, and in (|73p it depends on K, T, and the constant 5 from ()48p . Using the trace theorem 
[15] for the spaces H™ and the usual trace theorem for H m , from ([65]h (f68|) and (f73|) we get 

llfl , 5||ifi(sn T ) < C {lbi||i?2(an T ) + \\fi\x 1 =o\\H 1 (an T )} < C{||s , i||// 2 (dc T ) + [/i]2,*,t}- (74) 
||56[|ij-2(sn T ) < C { II (^3, 54) Hi^3(Of7 T ) + IK/ii) /i6)Ui=ollff 2 (9n T )} 

< C{l|(53,54)||i?3 (aQT) + ||/ii||h3 (Qt) . (75) 

Below we reduce problem (f54"J) - ([5T|) to that with (/n, /15, /i6, /17) = and (5,55,56) = 0, 
and in Sections [6] and [8] we will work with a reduced problem. If, however, we considered the 
original problem (f5"4"|) - (f5"T|) . taking into account (|72"j) - ([T5"j) . in the basic a priori estimate we 
should lose one derivative from f\ (in terms of H™ norms), one derivative from g\ = (51,52) 
and two derivatives from the source terms /n and 511 = (53,54) to the solutionis That is, the a 
priori estimate ()61[) is a roughened version of the basic a priori estimate. We prefer to sacrifice 
the lesser number of lost derivatives for obtaining such a reduced problem with homogeneous 
boundary conditions and /n = that for this problem we can derive an a priori estimate with 
no loss of derivatives from a new source term f\. In the future this will enable us to prove the 
existence of solutions of the reduced linearized problem by an almost classical argument (see 

mi). 

We first reduce problem ([54"]) - ([5?]) to that with fu = and 511 = 0. To this end we consider 
the auxiliary unknown V satisfying system (|55|) . the last two boundary conditions in (|56p with 
ip = and zero initial data: 

in 0, T , (76) 

on dn T , (77) 

(78) 

d 3 ip -d 2 (f> ^ 

-1 

1 


-x 
-x J 

5 The fact itself that we lose derivatives from the source terms to the solution is quite natural because the 
constant coefficients version of problem (|54jl - (|57p can satisfy the Kreiss-Lopatinski condition only in a weak 
sense. We prove this in Section [7] at least for particular cases. 



d t v - B^V + B 2 d 2 V + B 3 d 3 V = fn 



where 



E T2 - xU 3 = 53, E T3 + kH 2 = 54 
V = for t < 0, 

/ 



B> 



— >c 














—x 





-d 3 ip 








—x 


d 2 (p 





-d 3 (p 


d 2 (p 


—x 


d 3 ip 





1 





-d 2 (p 


-1 
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and 

(BiV, V) = U 3 E T2 - U 2 E T3 + U x (E T2 dw ~ E T3 d 2 <p) -x\V\ 2 . (79) 
Clearly, f) and e satisfy equations ([65]) : 

div~f) = /i 6! div~e = /i 7 in ft T - (80) 

Proposition 5.3 Problem ([75]) - ([75]) uni/i sufficiently smooth coefficients x, d 2 (f> and d^(p has a 
unique solution V G H s (Sl T ), with s > 0, /or a// (/ n ,flli) G H s+1 / 2 (n T ) x H s+1 / 2 (dn T ). This 
solution obeys the estimate 

\\V\\ H °(n T ) < c{l|/n||i^+V2(n T ) + [lanllir«+V2(an T )} , (81) 

where C > is a constants independent on the data (fn,gu) and depending on the coefficients^ 
and the time T. 

Proof. We note that the number of boundary conditions in (|77p coincides with the number of 
incoming characteristics of system (|76p . i.e., problem (|76j) — (|78j) has the property of maximality 
|llj . To prove the existence of solutions to problem ([76]) -(j78 p we reduce it to a problem with 
homogeneous and maximally dissipative boundary conditions. We first pass to a new unknown 
for which equations (|80p become homogeneous (with /ig = fu = 0). Actually, we even make 
system (|76p homogeneous. Let the vector- function V = (fi,E) which vanishes in the past 
satisfies system ([76]) and the boundary conditions 

E T2 = 0, E T3 =0 on dfl T - (82) 

Then, in view of {72]), (JS2J) and gSJ), we have 

(B 1 V,V)\ x i =0 = -h\V\ 2 >5\V\ 2 . 

That is, the boundary conditions (|82p are strictly (and maximally) dissipative and, hence, there 
exists a unique solution V obeying the estimate (see [TTl |2T] ) 

\\V\\h s (Q t ) + II^Ix^oIIh^^t) - C|l/nll-ff s (n T )- (83) 

The new unknown = V — V satisfies problem (|76|) -(|78 p with fjj = and g% and 54 
replaced by 

53=53 + ^3 and 54 = 54-^2- 
respectively. It follows from ([83]) that 

life, §4)\\h«(SIt) < C {ll/ii||H s (n T ) + lbii||H s (an T )} • (84) 
6 Below we use Proposition 15. 31 with an index s for which assumption (|43[) on the coefficients is enough. 
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Moreover, the new unknown V" satisfies the homogeneous equations f|S0() : 

div~ fjf = 0, div" e s = in Sl T . (85) 

Using the secondary symmetrization (|2ip of the Maxwell equations described in Section [21 
we now reduce the problem for V* to an equivalent problem with dissipative (but not strictly 
dissipative) boundary conditions. Taking into account (|2~T|). (f23j) and using ([85]) and the fact 
that 

i?idiv~ f) tt + i? 2 div" e tt = (lC 2 d 2 ip + /C 3 d 3 <£ - /Ci) d^V* + /C 2 <9 2 + /C 3 d 3 U tt , 
from system (|76p for V" (with /n = 0) we obtain 

BodtV* - BtdtV* + B 2 d 2 V i + B 3 d 3 V* = 



in fin 



(86) 



where £>i := B\{0) = B\ — x/?i — d 2 (pB 2 — d 3 (pB 3 . We choose 



Then 



Si 







V\ = 


X, l/ 2 


= ^3 = 0. 











—Kd 2 <$ 


—>cd 3 (p 





d 3 ip 


-d 2 (p 




—>cd 2 (p 


-2x 










-1-x 2 




—Kd 3 (p 





-2x 


d 2 (p 


1 + x 2 












d 2 (p 





—xd 2 (p 


—>cd 3 (p 




d 3 ip 
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Since |x| < 1 the hyperbolicity condition (f22j) holds, i.e., £>o > 0. It is noteworthy that from 
system ([86)) we can derive equations (j85j) and using them obtain the original system (jT6|h That 
is, systems ([76]) and ([86[) are equivalent and we may use ([86]) to prove the existence of a solution 

vK 

Assuming for a moment that g 3 = g^ = and using the homogeneous boundary conditions 
for 0, 



x^ 



J r 2 '"' l - 3 ; "r 3 

after straightforward calculations we get 



-xWy on 



(BiF»,F')U 1= o = 0. 

This means that the boundary conditions are dissipative (but not strictly dissipative) . To prove 
the existence of a solution V" we reduce the problem for it to that with homogeneous boundary 
conditions (with g 3 = g^ = 0). To this end we use the classical argument [T7] suggesting to 
subtract from the solution a more regular function satisfying the boundary conditions. Thus, 
referring to [T7j , we can conclude that there exists a unique solution V" , and this solution obeys 
the estimate 

IfT^fir) < C\\(g3,g4)\\ Hm +i/2( mT y (87) 
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Hence, there exists a unique solution V £ H s (Qt) to problem ([76]) — ([78]) and, in view of ([83]) . 
AMD and flHTJ), it obeys estimate □ 



Having in hand Proposition 15.31 in problem f|5-4|) — (|5T[) we can pass to the new "vacuum" 
unknown 

V = V - V. (88) 

Here to simplify the notations we again use the same letter for the new "vacuum" unknown as 
for the "vacuum" unknown in the linearized problem before passing to the "good unknown" V 
in (|5ip . The new unknown V satisfies the homogeneous system (|55p (with /n = 0) and the 
homogeneous last two boundary conditions in (|56p (with gn = 0). 

Concerning the "plasma" part of (|54 p - (j57p . the change of unknown (|88p only changes the 
source term in the second boundary condition in (|56p : 

q - (H, H) + (E, E) + [d iq ]<p = g 2 := g 2 + (H, U) - (E, E) on dQ T . 

By virtue of ([81 p and the trace theorem in H m , we have the estimate 

\\92\\H m {n T ) < C |||52||H m (9n T ) + ||Vja;i=ol|if m (an T )} 

< C {\\g 2 \\H™(dn T ) + ||^||h™+ 1 /2(q t )} 

< C{||5l|if™+i(an r ) + ll/n||if™+i(n T )} • ( 89 ) 

To reduce the "plasma" part of the linearized problem to that with homogeneous boundary 
conditions (with g\= gi = 0) and the homogeneous equations (f6T|) and (f62|) (with /15 = 55 = 0) 
we can exploit the same arguments as in for current-vortex sheets. Here we even describe a 
"shifting" function in an almost concrete form: 

t/ = (p,f 3 {i 1 ,0,0,i?,0), v x = X (x X )9x, p = X (x 1 )g 2 -^B 2 , 

5B 2 = i^(b,H) + (v,H)(H,u)-B 2 (v,u)^j , u = f 3 S, ^ = (^,0,0), 

where the cut-off function \ G Co°(^+) equals to 1 on [0, 1] and H solves the equation for H 
contained in (48) with v = v = (vi,0,0): 

d t H + (w,V)H-(h,V)v + Hdivw = f H + (h,V)v-HdiVi-v 1 d 1 H in Sl T , (90) 

where h = (H\ — H 2 d 2 (p — H^dsip, H 2 , H3) and fn = (/5,/6j/7)- It is very important that 
u>| x .i =0 = 0, i.e., the linear equation ([90]) does not need boundary conditions and we easily get 
the estimate 

[diH] s ^ : T + [H]s+i,*,t < C{[/i] s+ 2,*,r+ \\gi\\H"+ 2 (dn T )} ■ ( 91 ) 
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It is clear that U satisfies the first two boundary conditions in (|56p with g 2 = g 2 . Then the 
new unknown U = U — U satisfies the homogeneous boundary conditions ()56p and system (|54p 
with f\ = F, where 

F = (F 1 ,...,F 8 ) = fi-V' e (U,V+)U. 

In view of (p0|) . Fjj = (F5, Fq, Fj) = 0. Since the boundary conditions became homogeneous 
and Fjj = 0, it follows from Proposition 15.11 that the new "plasma" unknown U satisfies (|6ip 
and (ff)2"j) with /15 = 55 = 0. Moreover, taking into account (|59"P and (|9"Tj) . for the new source 
term F we have the estimate 

[F] S ,*, T < C[[/i] s ,*,t+ IK91, 52) ||if«+i(an T ) 

+ [<9i-H]s,*,t + [-ff]s+i,*,r| 

< C {[/l]«+2,*,T + ll/n||H s + 2 (n T ) + ll9llH s + 2 (an T )} • (92) 

In this paper we need ([92]) for s = 1: 

[F]i,*,t < C {[/i] 3 ,*,t + ||/ii||i/3 (nr) + \\g\\jjHdn T )} ■ (93) 

We write down the final form of our reduced linearized problem for the new unknown (U, V) 
and the interface perturbation <p: 



A d t U + J2AjdjU + CU = F in U T , 

d t V - B x d x V + B 2 d 2 V + B 3 d 3 V = in U T , 



d t <p = v N - v 2 d 2 ip - hd-np + ipdiv N , 

q = (H,H)-(E,E)-[d 1 q]ip, 

E T2 = d t (H 3 ip) - 8 2 {E W ) + kU^ 

E T3 = -d t {U 2 p) - d z {E w ) - x% 2 , 
(U, V, ip) = for t < 0, 



(94) 
(95) 



on 



(96) 



(97) 



where A a := A a (U) (a = 0,2,3), A\ := Ai(U,ip), C := C(U,<p). Recall also that solutions to 
problem (|94|) - (|97|) have the following very important property. Since Fjj = and equations 
and the boundary conditions (|96p are homogeneous, solutions to problem (|94p ~(|97 p satisfy 



div h = in fir, 

div~ f) = 0, div~ e = in CIt, 
Hn = H 2 d 2 tp + H 3 d 3 (p — (pdiHjy on dflx, 
Hn = 'H.ichp + T-Lsdzp — (p OiHn on 90t- 



(98) 

(99) 
(100) 
(101) 
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We recall that 

v N = vi - v 2 d 2 (p - v 3 d 3 ip, E Tj = Eidjtp + Ej, h = (H N , H 2 ,H 3 ), 

H N = Hi - H 2 d 2 - H 3 d 3 ip, U N =U\- U 2 d 2 (p - V. 3 d 3 ip, 
t) = (H N ,H2,H 3 ), e = {E N ,E 2 ,E 3 ), E N = E x - E 2 d 2 - E 3 d 3 0, 
and below we will also use the notations 

fi = (Hi,H Ti ,Ur a ), U Ti =Uxdi(p + U h <£ = (E 1; E T2 ,E T3 ). 

From now on we concentrate on the proof of the following basic a priori estimate for solutions 
to problem (f5ty -$57 $ : 

[U}l,*,T + \\V\\ HHciT) + M H HdU T ) < C[^]l,*,T, (102) 

where C = C(K,T) > is a constant independent of F. Taking into account (j93j) . (|8ip and a 
corresponding estimate for U, from f|102[) we obtain the desired a priori estimate (|60p for the 
original unknown (U, V) and the interface perturbation ip. It is worth noting that estimate (|102l) 
is with no loss of derivatives from the source term F to the solution. This will be useful for the 
proof of the existence of solutions to problem (fM|) - ([9T)) (this work is postponed to the future). 

6 Energy a priori estimate for the constant coefficients problem 

To clarify main ideas in a clearer form and to avoid unimportant technical difficulties we first 
prove estimate (|102p for the case of constant coefficients of (|94|) - (i97|) and dropped zero-order 
terms. In Section [8] we will extend this result to variable coefficients. That is, we now "freeze" 
coefficients of problem (|94p ~(|97 p and drop zero-order terms. Moreover, we assume that d 2 (p = 
d 3 ip. Then we obtain a constant coefficients linear problem which is, in fact, the result of the 
linearization of the original nonlinear free boundary value problem ([8]), (|13p . (|14p . (jl7p - (|19p 
about its exact constant solution 

U = U = (p,u, H, S) = const for x 1 > xt, 
V = V = (H, E) = const for x 1 < xt 

for the planar plasma- vacuum interface x 1 = xt, where x, with — 1 < x < 0, is a constant 
interface speed. We know that the components of this exact constant solution satisfy the relations 

(cf. g5D, flM}) 

0! = ux/f = m, H 1 =H 1 = 0, E 2 = xii 3 , E 3 = ->m 2 . (103) 
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Taking into account (|103p . we have the following constant coefficients problem: 

3 

A d t U + AjdjU = F in tt T , 

3=1 



(104) 
(105) 



(106) 



(107) 



d t V - B^V + B 2 d 2 V + B 3 d 3 V = in Sl T , 
d t ip = v\- v 2 d 2 (f - v 3 d 3 ip, 

q = n 2 (n 2 + *e 3 ) + u 3 {u 3 - xe 2 ) - EiE u 

E 2 = V. 3 d t f - Eid 2 ip + xH 3 , 

E 3 = -U 2 d t ip - Eid 3 ip - kK 2 , on 8Qt, 

(U, V,ip) = for t < 0, 

where A a := A a (U) {a = 0,2,3), A\ := A 1 (U) - >cAq, B x := B\ - xl, v = (u- (v,u)v) /f , etc. 
Moreover, solutions to problem (|104p - (|107p satisfy 

div H = 0, div" H = 0, div - E = in n T , (108) 

Hi = H 2 d 2 (p + H 3 d 3 if, %i = H 2 d 2 tf + H 3 d 3 (p on <9fi T . (109) 

From the physical point of view, the well-posedness of problem (|104p ~ (|107p can be interpreted 
as the stability of a planar relativistic plasma-vacuum interface (or the macroscopic stability of 
a corresponding nonplanar interface). 

For problem fllQ4f) — ( TlOT]) we now derive the a priori estimate (I102p under the sufficient 
stability condition that the constant fi = \fi\, with 

A = E% + V 2 H 3 - t>3%2, 

is small enough (cf. (|59p ) and the vectors (H 2 , H 3 ) and ('H 2 ,T-L 3 ) are nonzero and nonparallel to 
each other (cf. (i58j) ). i.e., 

H 2 H 3 - H 3 H 2 + 0. (110) 

Using the secondary symmetrization ([2T]) of the Maxwell equations, as in ([86]) . we equiva- 
lently rewrite system (11051) as 



B Q d t V-BidiV + B 2 d 2 V + B 3 d 3 V = in U T , 
where Bi = Bi — xBq. Now our choice of Ui is the following: 

y = V = (t!i,V 2 ,V 3 ) 



(111) 



(112) 



(vi = x, see (|103p ). Since \v\ < 1 the hyperbolicity condition ()22[) holds. By standard arguments 
of the energy method applied to the symmetric hyperbolic systems (11040 (multiplied by 1/r) 
and (jllip we obtain 

m + 2 Ln Qdx ' ds ^ C (\\ F \\Un T ) + J*Hs)a s y (113) 



2G 



where 

I(t)= [ (l(A U,U) + (B V,V))dx, Q = hB 1 V,V)\ xl=0 -^(A 1 U,U)\ xl=0 . 

JrI v r ' 2 2r 

Taking into account (|35|) . we write down the quadratic form Q: 

Q = {-k(H\ +%l + El + El) + Ui(v 2 U 2 + v 3 U 3 ) + E x (v 2 E 2 + V3E3) 
+ hHx(v 2 E 3 - v 3 E 2 ) + (1 + x 2 ){U 3 E 2 - H 2 E 3 ) - qvi}\ xl=Q . 

Using the boundary conditions (|106p and the second condition in (|109|) . after relatively long 
calculations, which are omitted here, we get 



Q = pL {Exdttp + (n 2 + xE 3 )d 3 ip - (H 3 - xE 2 )d 2 ip}\ xl 



x 1 ^ 



or 



Q = dt (A^iL.i=o) + d2 (w(xE 2 - U 3 )\ x i =0 ) + d 3 (MH 2 + xE 3 )\ xl=0 ) + Qo, 
and in view of the fourth equation in ()106|) and the third equation in ()109|) , 

Qo = {dtEi + xd 1 E 1 + d 3 U 2 - d 2 H 3 )\ x i =0 = 0- 
Then, it follows from (|1 13j) that 

I(t) + 2 J fOpExl^dx' < C (\\F\\1 2{Qt) + J I(s) ds) . (114) 

At first sight inequality f j 1 1 4 j) is absolutely useless. Indeed, we cannot apply it for deriving an 
L 2 a priori estimate. On the other hand, if we differentiate systems (|104|) and (jllip with respect 
to x 2 , x 3 and t, we easily obtain the following counterparts of (|114p for the first-order tangential 
derivatives of U and V: 

I a (t) + 2 J^fidavdaE^i^dx' < c(\F%<r + j I a (s)dsY (115) 

where a = 0, 2, 3, 

Ut)= [ (UA d a U,d a U) + (B d a V,d a V))dx and d := d t . 

The terms d a i~pd a E\ \ x x = q appearing in the boundary integral are actually lower-order terms and 
below we explain how to treat them by using our important assumption (jllOp and passing to 
the volume integral. 

Thanks to (fTTOl) . from (PT091) we get 

r, (H3H1 - H 3 Ui)\ x i =0 1 x 
2 tp - -^—^ - a^i^i^ + a 2 Hi\ x i =0 , 

H 2 n 3 -H 3 n 2 (n6) 

[% 2 H 1 - H 2 %i)\ x i = q 2 2 

C 3 </? = ^-^ = a^ila-i^ + a 2 Hi\ x i =0 . 

ri 2 fi 3 — ri 3 fi 2 
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First we consider spatial tangential derivatives in (|115p . i.e., let a = k with k = 2 or k = 3. By 
using (|116p one gets 



<9^d fc £iU =0 = laZHtdkEi + aZnidkEi) 1 , k = 2,3. (117) 

Concerning the term dt<p dtEi\ x i =0 , it can be reduced to the multiplication of spatial tan- 
gential derivatives of U and V. Indeed, (|116p and the first boundary condition in f)106[) imply 

d t ip = a\H x \ x i = Q + a%Hi\ x x =0 + vx\ x x =G , (118) 

where the constants a® can be easily written down. It follows from the fourth equation in (|106|) 
and the third equation in (|109p that 

d t E x = d 2 {U 3 + *E 2 ) - d 3 {U 2 - xE 3 ). (119) 

Then, (fTT7V(fTT9|) yield 

3 

(dtpdtE! + d 2 pd 2 E 1 + d 3 tpd 3 Ex)\ x i =0 = 5Z(M fc y,a fc y)Ui =0 , (120) 

k=2 

where Y = (vx, Hx,Hx) and are rectangular matrices which constant coefficients depending 
on (U, V) are of no interest. 

To treat the integrals of the lower-order terms like H\dkE\\ x i = Q contained in the right-hand 
side of (|120p we use the same standard arguments as in [251 [261 12Z]- That is, we pass to the 
volume integral and integrate by parts: 

f Hxd k Ex\ xl=0 dx' = - f d 1 (H 1 d k E 1 )dx= [ {d h HxdxEx-dxH x d k Ex}dx. 

JR 2 JR3_ JR3_ 

After that the normal (x 1 -) derivatives of vx, Hx and V can be expressed through tangential 
derivatives by using the facts that the boundary x 1 = is noncharacteristic for the "vacuum" 
system (|105j) and v\ and Hx are noncharacteristic "plasma" unknowns. Namely, since det-Bi 7^ 0, 

dxV = (Bx)-\d t V + B 2 d 2 V + B 3 d 3 V), (121) 

and taking into account ([35[) . (|37p and the first equation in (|108j) . we have 

dxvx \ 1 f t / — - - \ 1 

= \J T [A d t U + A 2 d 2 U + A 3 d 3 U + F)\ , (122) 
[1/T)dxq J \ /J (12) 

OxHx = -d 2 H 2 - 8 3 H 3 , (123) 
where {• • • }( 12 ) denotes the first two components of the vector inside braces. Hence, we obtain 
2/t V / d a ipd a Ex\ x i =0 dx' = fi / (QZ,Z)dx + {QxF,Z)dx, (124) 
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where Q is a quadratic matrix of order 42 which elements can be explicitly written down if 
necessary, Z = (dfU, dfV, d2U, 82V, d^U, 83V) , and Qi is a rectangular matrix which elements 
are of no interest. 

Now, using (|124p . we consider the sum of inequalities ()1 15[) : 

f ((a + £Q)Z,Z))ds + ! (Q 1 F,Z)dx<C I [Flla + [ Ia(s)ds\, (125) 

where 21 = diag ((1/r )Aq,Bq, ... ,(1/T )Aq, Bq) is the block-diagonal positive definite matrix of 
order 42. If n = \ is small enough, then the matrix 21 + fiQ is also positive definite (this is 
enough to derive estimate ([102p . see below). In principle, we can analyze the sufficient stability 
condition 

2l + /iQ>0 (126) 

if necessary. Since the matrix 21 + fiQ is of order 42 and in RMHD the matrix Aq is rather 
complicated, this is almost impossible to be done analytically, but a corresponding numerical 
analysis of condition (|126p seems straightforward. 

In this paper we restrict ourself to the assumption that fi is small enough. Using the Young 
inequality and the elementary inequality 

\\nt)\\l n) <C[F]l tt (127) 

following from the trivial relation 

i||F( t )lll l(R , ) =2/^(F, 8( F)dx, 

we estimate: 



7k- 



(QiF, Z)dx <c{e \\Zml^ +) + \ \\F\\1*,t} , (128) 



where e is a small positive constant. Choosing e to be small enough and using ()126p and (|128p . 
from inequality (|125p we obtain 



WZWWUpV ^ C {[F]irt + \\Z\\t 2{nt) ) . (129) 

Here we can already use Gronwall's lemma, but we prefer to do not do this and follow arguments 
which can be then extended to the case of variable coefficients (see Section [8]). 

We can easily include the L2 norms of U, V, ad±U and d±V in inequality (|129p . Indeed, to 
estimate the weighted derivatives of U we do not need boundary conditions because the weight 
cr\ x i =0 = 0. By applying to system (|104p the operator ad\ and using standard arguments of the 
energy method, we obtain the inequality 

hdiUit)\\l 2 ^? +) <C{[F\i*? + P\i*,t} (130) 
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(actually for the present case of constant coefficients the norm [U]\^ t in (|130p can be even 
replaced by ll^l^ll^m^)- It follows from (|12ip that 

\\diV{ml 2 ^ +) <C\\V\\ 2 m{nt) . (131) 

Summing up inequalities ()129j) — f)131|) and the elementary inequalities for U and V like ()127p . 

\\U(t)\\l 2 ^ + ) < Cp?w \\V(t)\\l H) < C\\V\\ 2 m{nt) , (132) 

we get 

+ \\\V(t)\\\ 2 HHRl) < C (Win + [U]l^ t + \\V\\ 2 HHnt) ) , (133) 

with 

W\V(t)\\\ 2 HHRl) ■= \\V(t)\\mo* +) + \\dtV(t)\\l 2{RS+) . 
It remains to estimate the interface perturbation (p. By using f|122|) and (|123p we first 
estimate the noncharacteristic part U n = (q,vi,Hi) of the unknown U: 

II^(*)IIL (R 3) < C(K){[F]\^ T + \\\U{t)\\\i*} ■ (134) 
Then thanks to the trace theorems, (|13ip and (|134p imply 

\\UnU=o{t)\\l 2{ ^) < C{[F}\^ T + |||{/(i)|||t} , (135) 

rU=o(i)llL( K2) < C\\V\\ (136) 

The multiplication of the first boundary condition in (|106|) by 2tp, the integration over the 
domain M 2 , and the usage of (|116p . (|118p . estimates (j!35|) and (j!36|) result in 

||b(t)||||i {K2) < C{[F]\^ T + [U\l #tt + \\V\\ 2 Hl(Qt) + y\\l 2(dnt) }. (137) 
Summing up inequalities (j!33|) and ()137|) we obtain 

III^WIII?,* + \\\vm\l H ml) + Mt)\\\ 2 H ^) 

<c{[F} 2 ^ T + [U} 2 ^ t + \\V\\ }• (138) 

Applying Gronwall's lemma, from (|138|) we derive the desired a priori estimate (|102|) . 

Remark 6.1 In fact, for problem ()104p ~ (jl07p we have constructed a dissipative 1-symmetrizer 
|25j . The index 1 means that we prolong our hyperbolic system up to first-order derivatives 
and manage to construct a secondary symmetrization of the prolonged system with dissipa- 
tive boundary conditions (see |25j). The divergence constrains (|1Q8|) play a crucial role in this 
construction. Note that in |24| we constructed a dissipative O-symmetrizer and proved an L2 es- 
timate for the constant coefficients linearized problem for current- vortex sheets and then derived 
an estimate in Hi for the variable coefficients problem. Here, for the plasma-vacuum interface 
problem we have to prolong the hyperbolic system up to first-order derivatives already on the 
level of constant coefficients. 
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7 Normal modes analysis for particular cases 



There appears the following natural question. Can problem (|104|) — (jlOT[) be really ill-posed or 
just our approach in Section [6] does not allow us to derive an energy a priori estimate for the 
whole range of admissible parameters of this problem? Below we show that this problem can 
be indeed ill-posed. In other words, unlike planar fluid-vacuum boundaries [12\ |2"7] and planar 
nonrelativistic plasma- vacuum interfaces [28J, a planar relativistic plasma- vacuum interface can 
be violently unstable, i.e., Kelvin- Helmholz instability can occur. Mathematically this means 
that there is a range of admissible parameters of problem (|104p fjlQ7|) for which the Kreiss- 
Lopatinski condition [TOl Q3] is violated. 

Unfortunately, due to principal technical difficulties a complete Fourier-Laplace analysis of 
problem ()104|) — (j!06|) seems impossible (even numerically). Therefore, now our goal is just to 
find particular examples of the ill-posedness of problem (|104|) — f|106[) . We restrict ourself to a 
case when the structure of the matrices A a in ([3D is not so complicated and we can perform 
analytical calculations. This is the case when the 3-vectors of velocity and magnetic field are 
collinear and have only one nonzero component. Such a case was already considered in |23j for 
relativistic MHD shock waves. 

The test of the Kreiss-Lopatinski condition is equivalent to the usual normal modes analysis, 
i.e., the construction of an Hadamard-type ill-posedness example. We seek exponential solutions 
to problem f|104[) (|106[) with F = and special initial data: 

U = Uexp{Tt + t p x 1 +i(j',x')} , (139) 
V = V enp {rt + ^x 1 +i(V,x')} > ( 14 0) 
(p = <pexp{Tt + i('y , ,x')}, (141) 

where U and V are complex-valued constant vectors, r, £ p , £ v and (p are complex constants, and 
Y = (72)73) with real constants 72,3. The existence of exponentially growing solutions in form 

TOHBD, with 

5ftr>0, 3?£ p <0, K£ v <0, (142) 
implies the ill-posedness of problem (|104p - (jl06p because the consequence of solutions 

(U(nt, nx), V(nt, nx), (p(nt, nx')) exp(— \/n), n = 1,2,3,... , 

with U, V and tp defined in (I139p - (|14ip . is the Hadamard-type ill-posedness example. Since 
the last equation in f| 104[) for the entropy perturbation S plays no role in the construction of an 
ill-posedness example, we will suppose that U = (p,u,H). 

Theoretically, we can construct a ID ill-posedness example (with 7' = 0) on a codimension-1 
set in the parameter domain. But, this is not the case for problem (|104p - (|106p . Indeed, in ID 
we have E\ = and the energy inequality (jll4j) . which can be obtained even for x > 0, implies 
well-posedness. 
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Thus, we may assume that 7' 7^ and even, without loss of generality, |7'| = 1. Moreover, in 
view of the big complexity of the RMHD system, we will construct a 2D ill-posedness example 
with 73 = 0, i.e., 72 = 1. Let also consider the particular case v% = i>2 = and H2 = for the 
relativistic magnetoacoustic system (|114|) (recall that H\ = 0, see (|103p ). Since v\ = x, we have 
x = 0. However, by continuity the below ill-posedness examples can be extended to the case 
\x\ <C 1. For the particular case v\ = Vi = Hi = H2 = the matrices A a (without the last rows 
and columns) take the following "visible" form: 

/ f/(pa 2 ) 

phf + HVf 

phf + Hl/Y 
A n 



\ 

/o 
1 










ij 3 




1 









phT + Hl/T 







Hi 




















H* 



An 



V3 




ph/f 





/ 



1 


















































1/f 













1/f 













f 


) 



J 



1 











^3 



V 




















/ 



H 3 

0). Since 73 = 0, the matrix A3 is of 



(see also [23 for the analogous case vi = V3 = H2 
no interest and we do not write down it here. 

Substituting ([139]) and (jTiOl) into (fTQlD (with F = 0) and ([105]) respectively, we get the 
dispersion relations 

det(rA + £ p li + iA 2 ) = 0, det(r/ - £ v #i + iB 2 ) = 
from which, omitting technical calculations, we find the unique roots 

Vl + rr 2 , 



with property (I142D . where 

w 2 (f 2 + m 2 c 2 s ) 



Vl + r 2 , 



(143) 



ur 



1 - vie 2 > 



H3 



m 



c 2 {l + m?w 2 ) ' " ~ " "" a^p 

Consider first the particular case H2 = 0. For this case our basic assumption (jllOp is violated. 
Substituting (|139)) — (|141|) into the boundary conditions (|106|) . we obtain 



Tlf = Vl, 



H3H3 - EiEx, E 2 = tH 3 <p - iEup, 
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and excluding the constant <p we come to the relations 

q + EtEx - n 3 n 3 = 0, (iE x - rn 3 )vi + tE 2 
Prom (fT39|) and (fT10|) we have 



tpQ + T\phT + 
Then, (fTiij) and (fT13|) imply 



0, rE 1 =iH 3 , iE 2 = ^E x . 



(144) 



(145) 





1 





Ei + irHa 


\ 


( q \ 




Cp 


r(phf + HI/T) 







Vl 


V 





iE x ~ tH 3 




) 


\Ei ) 



(146) 



Since we are interesting in solutions with (U, V) ^ 0, the determinant of the matrix in (|146|) 
should be equal to zero. Taking into account (|143|) . this gives us the following final equation for 
t: 



rVl + r 2 = 1{E X + i^3T)Vl + rr 2 , 



(147) 



where 



phT 2 + Hi 



>0, 



r > 1. 



If Ei = Hz = 0, then /i = and equation (|147p has only roots with SRr = 0. That is, at least 
for the 2D problem (recall that 73 = 0) the Kreiss-Lopatinski is satisfied in a weak sense. Note 
that for the case /i = we can derive an L 2 a priori estimate for problem (|lU4|) - (|107p (in 3D) 
even if assumption (jllOp is violated. This follows from (|115p . 

If "%3 = and Ei ^ 0, then graphically it is easy to see that equation (|147p has a root r 2 > 0, 
i.e., we have an "unstable" root r > 0. Clearly, there is an "unstable" root r with Sir > at 
least if 1 7^3 1 is small enough. It seems that one can prove that (|147p has a root with property 
(|142p even for arbitrary values of I-L3. We do not study this technical problem here. We only 
note that for the case r — > 1 which correspond to the ultrarelativistic limit c s — )■ 1 equation 
(|147p has a root with 3?r > for arbitrary %3 and E\ 7^ 0. Thus, we see that if assumption 
(jllOp is violated, then the linearized problem can be ill-posed. 

We now consider a particular case when assumption (jllOp holds. Let 7i 3 = and H-yH. 2 7^ 0. 
Substituting (|139|) - (|141|) into the boundary conditions (jlUtjp . after the exclusion of constant (p 
we obtain 

q + E1E1 - U 2 U 2 = 0, iEiVi + tE 2 = 0, t7 2 wi + E 3 = 0. (148) 
It follows from (fTQ4"P and ^05]) that 



( p f+r ^f + ^ D 1= 0, £ v E 3 = Tn 2 , iE 2 =£ v E 1 . 



(149) 
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The algebraic system (|148p . (|149p has a nonzero solution (q,vi,E) if 



T 



Vl + r 2 = - U\{\ + r 2 )) Vl + rr 2 



(150) 



(we omit intermediate simple calculations) or G(y) = 0, where y = r 



2 and 



= yVlT^ - £{Ef - n\{\ + y)) ^/T+7y~. 
If > Hi, then G(0) < and G(y*) > 0, where y* = (Ef/Hj) - 1 > (recall that we 



then equation (|150p has an "unstable" root r > 0, i.e., problem (|104p (jlQ7j) is ill-posed. If 



at least for r — > 1 it is easy to see that equation (|150p has only roots with 9?t = 0. That is, 
again the Kreiss-Lopatinski can be satisfied only in a weak sense. 

For the particular case under consideration, n = \E\ — v^H.2\- In Section [6] we could not 
unfortunately concretize how small should be fi to guarantee the fulfillment of (|126|) . Clearly, 
if Ei is close to v^hlii then (|126p holds. If, for example, = 1/2, we see that the condition 
Ei ~ ^3^2 is more restrictive than (|152p . In general, if Ei is close to V3I-I2, the instability 
condition (|151 1) is violated. 

8 Variable coefficients analysis 

The principal difference of the derivation of estimate (|102p for the variable coefficients case in 
comparison with the constant coefficients one is the appearance of additional "not dangerous" 
lower-order terms and the fact that the usage of spaces H™ is a principal pointQ We now give 
short comments about the adaptation of arguments from Section [6] to the variable coefficients 
problem (p4"l) - (l97jl . 

First of all, for problem AMD"© we can easily get estimates (fT30]) - (fT32|) . (fT34l - (fT36|) . where 
now U n = (q, vn, Hn)- For example, the usage of the norms of Hi in the right-hand side of 
(|134|) is really necessary because the boundary matrix J T AiJ = 612 + M (see ([55]) . (JHTJ) ) and 
we can guarantee that the matrix N is zero only on the boundary: N\ x i = 0. Concerning the 
preparatory estimate (|137p . by virtue (|100p . ([10ip and assumption ([58j) . we can also obtain it 
because now in (|116p . (|118p we just have additional zero-order terms in ip and should replace 
Hi by H N , %i by U Nl and v x by v N . 

7 For the constant coefficients problem we could estimate the "plasma" unknown U in the norm ||?7|| H i = 



consider the case H2 7^ 0). Hence, if 




(151) 




(152) 



\U\\ L v + \\dtU\\ L2 + \\ckU\\ L2 + \\d 3 U\\ L2 . 
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Using {[99]) , we get system ([86]) for V and our choice is again (|112p (but now Vj are not 
constants). Then, we deduce the following energy inequality for d a U and d a V with a = 0, 2, 3: 

/«(*) + 2 / Q a dx'd S <c{[F]^ T + [[/]2 +||y||^ }, (153) 
where I Q were determined just after (|115p and 
<2« = l(Bid a V,d a V)\ x i^o - i(li9 a C/,5 a f/) 1^=0 

= ~(Bld a V, d a V)\ x l =0 - daqdaVNlx^O- 

We first calculate the quadratic form Q determined in Section [6] (see (|113p ): 

Q = {v 2 UiH 2 + vzHiHz + {v 2 d 2 (p + hdz<i>)Wi + qv N }\ x i =Q ■ 

This is a point when we have to perform big technical calculations. Omitting these rather long 
calculations which are really straightforward, we obtain 

Q = A {E N d t <p + (H T2 + xE 3 )d3(p - {H T3 - xE 2 )d 2 <p}\ x i =0 + l.o.t., 

where l.o.t. is the sum of lower-order terms like [diq] vn\ x i = q V- While obtaining the above 
formula for Q we used assumptions ([45]) . the boundary conditions ([96]) and (llOip . Taking into 
account the equation 

d t t - V~ x Sj + xdie + B(<p)t = 
following from ([95]) and applying the second equation in ([99]) . we rewrite Q as follows: 

Q = d t (fupE N \ x i =0 ) + d 2 (p,<p(xE 2 - n T3 )\ x i =0 ) + d 3 {fnp{Ur 2 + xE 3 )\ x i= ) + l.o.t. 

Then we have 

/ Q a dx'ds= / fid a ipd a E N \ x i =0 dx' + / (l.o.t.) da/ds, (154) 
Jdfit Jm: 2 JdQt 

where l.o.t. is the sum of lower-order terms like 

cd a v N \ x i =0 d a (p, cd a q\ x i =0 d a (p, cd a Hj\ x i =0 d a (f, cd a Ej\ x i =0 d a ip, 

etc., here and below c is the common notation of a coefficient depending on the basic state (j42p . 
By using the variable coefficients counterparts of (|116|) and (|118|) (where we have additional 
zero-order terms in (p; H\ is replaced by Hn, etc.), we reduce the lower-order terms above to 
those like 

cHjsrd a visr\ x i =0 , cT-L^daHj^^, cH]sfd a Ej\ x i =0 , cHNd a q\ x i=o, 

etc., or even "better" terms like c(pd a Ej\ x i =0 . By integration by parts such "better" terms can 
be reduced to the above ones and terms like c(pEj\ x i =0 . Actually, if a = 0, then the integration 
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by parts leads to the appearance of terms in the form dt(- ■ ■ ) which do not disappear after the 
integration over the domain d£lf hi this case we use the Young inequality (see explanations 
below). 

We treat the terms like cHNd a VN\ x i =0 by passing to the volume integral and integrating 
by parts (the only difference from the analogous manipulation in Section [U] is the appearance of 
variable coefficients): 



Jn t 

where c| x i =0 = c. If a = 2 or a = 3 the last integral above is equal to zero. But for a = we 
have: 



Using the Young inequality, (|132p and (|134p (with U n = (q, vn, Hn)), we estimate the last 
integral as follows: 



where e is a small positive constant. Concerning the lower-order terms ccpEj\ x i =0 , cipHj^\ x i = Q, 
etc., we treat them by using the trace theorems in H 1 and H}. 

If q = 2 or a = 3, we apply to the boundary integral in (1154p the same arguments as for 
the last integral there and take into account (|13(J|) - (|132p . (|134p - (|137p . If a = 0, we first use the 
variable coefficients counterparts of (|118p . (|119p and then apply the same arguments as for the 
cases a = 2 and a = 3. In particular, for variable coefficients in (|119p we have En instead of 
Ei, H T2 instead of H2, and H T3 instead of ^3. 

Thus, choosing e and the constant e from the variable coefficients counterpart of (|128p to be 
small enough, we finally get inequality (|138|) . provided that the constant fi* from Theorem 14.11 
is also small enough^! Then, applying Gronwall's lemma, we derive the desired a priori estimate 
(|102p . We did not discuss assumption (|43p . but we can carefully check that this assumption 
on the coefficients is indeed optimal. In this connection, we refer to [21] where the analogous 
assumption appeared for the linearized problem for current-vortex sheets. 

8 Actually, for variable coefficients the less restrictive sufficient stability condition is the same as in (| 126[) . The 
only difference is that coefficients of 21 and Q axe now functions depending on U and V. 





cd a H^diVN + (d a c)H^fdiVN — cdiH^d a VN — {di^H^diVN fdxds 
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9 Concluding remarks 



We have derived a basic a priori estimate for the linearized plasma-vacuum interface problem 
in RMHD under a sufficient stability condition. This condition precludes the violent instability 
of planar relativistic plasma- vacuum interfaces. We have also proved that, unlike the situation 
for fluid- vacuum free boundaries \1'2\ |2"7] and nonrelativistic plasma- vacuum interfaces [28] , the 
constant coefficients linearized plasma- vacuum problem in RMHD can be ill-posed due to a large 
enough vacuum electric field, i.e., Kelvin-Helmholz instability can really occur. 

In the basic a priori estimate (|60p obtained for the variable coefficients linearized problem we 
have a loss of derivatives from the source terms to the solution. Therefore, we plan to prove the 
existence of solutions to the original nonlinear problem by a suitable Nash-Moser-type iteration 
scheme. We suppose that the scheme of the proof will be really similar to that in [61 [26l 127] . 
The important point in this direction is the derivation of a so-called tame estimate for solutions 
(U,V,ip) of the linearized problem in H^(Qt) x H s (£It) x H s (8Qt) for s large enough. 

In this paper we restricted ourselves to the case of special relativity. At first sight, to cover 
general relativity we could try to follow a scheme proposed in [27J for the fluid-vacuum problem. 
Namely, in [27] we proposed to resolve the matter equations by Nash-Moser iterations whereas 
at each Nash-Moser iteration step we find the metric from the Einstein equations by Picard 
iterations. It is quite possible that such a scheme of the proof of the existence of solutions to the 
nonlinear problem could work (however, we are not yet absolutely sure in success even for the 
general relativistic fluid- vacuum problem) , but it is still unclear how to prove the uniqueness of 
a solution. 

The point is that to prove uniqueness we have to derive an priori estimate for the linearized 
problem with all the lower-order terms. It means that we have to pass to the Alinhac's "good 
unknown" even in the linearized Einstein equations written in the form of a symmetric hyperbolic 
system [18]. That is the linear boundary conditions [d 7 g a p] = for the metric g (see [27]) after 
the passage to the "good unknown" g a p become 

[d y g a p] = -<p[did^g a p], 

where g a p is an "unperturbed" metric and [•] denotes a jump at x 1 = 0. The right-hand sides in 
the above boundary conditions can be nonzero and the fact that for the fluid-vacuum problem 
the symbol associated to the free surface is not elliptic prevents obtaining the a priori estimate. 

It is astonishing that for the general relativistic plasma-vacuum interface problem we do 
not have such a difficulty at all because there is a natural stabilizing physical mechanism. This 
mechanism is the magnetic field (!). Indeed, under assumption (|58p the interface symbol is elliptic 
and, therefore, the lower-order terms appearing in the boundary conditions for the metric g after 
the passage to the "good unknown" are not dangerous anymore. Thus, in the future after the 
proof of a nonlinear existence theorem for problem ([5]), (fTU|) . (fT3|) . ([BP , CTJ, (f!8j) we hope to 
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extend this result to the general relativistic version of this problem. 
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